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Abstract.  A  frequency-domain  distributed  temperature/strain  sensor  based  on  a  longitudinally  graded  optical 
fiber  (LGF)  is  proposed  and  evaluated.  In  an  LGF,  the  Briilouin  scattering  frequency,  vB,  changes  (i.e.,  is  chirped) 
lengthwise  monotonically  and  thus  every  position  along  the  fiber  has  a  unique  vB.  Any  change  locally  (at  some 
position)  in  the  fiber  environment  will  result  in  a  measurable  change  in  the  shape  of  the  Briilouin  gain  spectrum 
(BGS)  near  the  frequency  component  mapped  to  that  position.  This  is  demonstrated  via  measurements  and 
modeling  for  an  LGF  with  local  heating.  The  LGF  is  one  with  ~100  MFIz  Briilouin  frequency  gradient  over 
16.7  m,  with  1.1  and  1 .7  m  segments  heated  up  to  40  K  above  ambient.  A  measurement  of  the  BGS  can  enable 
the  determination  of  a  thermal  (or  strain)  distribution  along  a  sensor  fiber,  thus  rendering  the  system  one  that  is  in 
the  frequency  domain.  A  sensitivity  analysis  is  also  presented  for  both  coherent  and  pump-probe  BGS  meas¬ 
urement  schemes.  The  modeling  results  suggest  that  the  frequency-domain  systems  based  on  fibers  with  a 
chirped  Briilouin  frequency  are  highly  suited  as  inexpensive  event  sensors  (alarms)  and  have  the  potential 
to  reach  submeter  position  determination  with  sub-1 -K  temperature  accuracies  at  >1  kHz  sampling  rates. 
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1  Introduction 

The  Briilouin  scattering  frequency,  vB,  is  related  to  the  opti¬ 
cal  wavelength,  A0,  and  the  acoustic  velocity,  V,  by 
vB  =  2Vnmode//l0.  It  is  well  known  that  both  the  acoustic 
velocity  and  the  refractive  index  of  a  glass  usually  are  func¬ 
tions  of  its  thermomechanical  environment;  i.e.,  temperature 
and  any  applied  stress  or  strain  to  the  material.  The  word 
“usually”  is  used  as  a  qualifier  because  it  is  possible  to  real¬ 
ize  multicomponent  materials  of  which  physical  properties, 
such  as  the  refractive  index,  are  immune  to  temperature  or 
strain/stress. 1-4  However,  in  conventional  pure  silica  or 
GeCA-doped  silica  core  optical  fibers,  the  Briilouin  scatter¬ 
ing  frequency  is  dependent  on  both  the  temperature  and 
strain,  with  the  usual  dependence  being  an  increasing  fre¬ 
quency  with  increasing  temperature  or  strain.6 

Since  the  temperature  and  strain  dependencies  of  the 
Briilouin  frequency  (typically  ~1.2  MHz/K6  and 
500  MHz/%,7  respectively,  near  1550  nm  for  conventional 
fiber7)  are  not  insignificant,  it  is  natural  that  this  phenomenon 
has  been  adapted  for  use  in  distributed  sensor  technologies,8"1 1 
with  very  impressive  results.  High-sensitivity  and  resolution 
measurement  schemes  and  configurations  have  been  demon¬ 
strated,  1 2"2"  and  a  commercial  market  exists  for  such  distributed 
sensing  systems.  Resolution  is  defined  in  this  context 
to  be  the  ability  to  resolve  two  events  that  are  in  close  proximity 
to  each  other,  or  alternatively,  to  resolve  a  certain  temperature  or 
strain  change  within  a  given  range  in  the  sensor  fiber. 
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Sensitivity  is,  then,  the  minimum  detectable  change  in 
temperature  or  strain.  In  recent  years,  there  has  been  a  signifi¬ 
cant  push  toward  the  achievement  of  remarkable  spatial  and 
thermal/strain  resolutions  (reviews  may  be  found  in  Ref.  21 
and  22).  However,  in  many  applications,  sub-centimeter  meas¬ 
urement  resolutions  (or  rather  spatial  resolutions  of  order 
«;10_4  of  the  total  sensor  length)  are  not  necessary,  and  instead 
robustness,  simplicity,  fast  sampling  times,  and  low  cost  are 
primary  concerns.  Examples  of  such  cases  might  be  fire  detec¬ 
tion  in  buildings,  tunnels,  etc.,2,  leak  detection,24  machines,25 
and  perimeter  security  (or  intrusion  sensing),26  to  name  a  few. 

Such  Briilouin  systems  are  labeled  as  “distributed” 
because  they  are  capable  of  correlating  a  lengthwise  thermal 
or  strain  distribution  in  a  deployed  sensor  fiber.  In  a  way 
fully  analogous  to  a  light  detection  and  ranging  (LIDAR) 
system,  time-domain  measurements  are  achieved  by  launch¬ 
ing  a  pulse  of  narrow  linewidth  light  into  a  fiber  and  receiv¬ 
ing  time-gated  returns  corresponding  to  Briilouin  scattering 
or  gain  from  known  positions  along  the  fiber.  Thus,  these 
systems  are  commonly  referred  to  as  Briilouin  optical 
time-domain  reflectometry  or  analyzer  systems,14'27  respec¬ 
tively.  Some  work  has  been  done  in  the  frequency  domain 
(i.e.,  Briilouin  frequency-domain  systems28,29  as  well  as 
non-Brillouin  configurations30),  with  promising  and  impres¬ 
sive  results,  but  Briilouin  frequency-domain  systems  are  far 
rarer  than  their  time-domain  counterparts. 

Continuing  with  the  LIDAR  analogy,  typical  atmospheric 
systems  are  sometimes  employed  to  measure  weather 
conditions,  such  as  temperature  or  wind  speeds.  '1  They  rely 
on  atmospheric  scattering  phenomena,  such  as  Rayleigh 
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scattering32  or  resonance  fluorescence,33  as  a  probe  from 
which  atmospheric  attributes  may  be  determined.  With  a 
LIDAR  system,  one  is  limited  to  designing  optical  systems 
based  on  the  whims  of  Mother  Nature.  One  would  be  hard- 
pressed  to  design  and  implement  a  custom  atmospheric  col¬ 
umn  that  would  lead  to  a  more  sensitive  probe  of  the 
environment. 

While  many  systems  utilize  conventional  optical  fiber  in 
distributed  sensing  applications,  this  is  tantamount  to  sub¬ 
mitting  to  Mother  Nature’s  whim,  and  there  is  nothing  nec¬ 
essarily  precluding  the  design  of  a  more  ideal  optical  fiber 
that  can  serve  to  improve  the  sensitivity  or  to  enhance  the 
capabilities  of  a  Brillouin-based  distributed  sensor.  For 
example,  acoustic  waveguide  designs  can  be  implemented 
that  give  rise  to  multiple  Brillouin  frequencies  being  present 
in  a  single  fiber  such  that  simultaneous  measurements  of 
temperature  or  strain  can  be  made34  ’3  with  a  single  sensor 
fiber.  More  recently,  it  was  proposed  that  designer  materials 
can  be  engineered  for  which  Brillouin  frequencies  are 
immune  to  temperature’8  or  strain.  ’7  One  may  even  envision 
the  possibility  of  a  material  with  enhanced  Brillouin  gain. 

Also  recently,  a  novel  way  of  fabricating  a  fiber,  for  which 
material  composition  is  lengthwise  graded  (longitudinally 
graded  fiber,  or  LGF),  was  proposed  and  experimentally 
demonstrated.’8  Briefly,  a  core  slug  is  drilled  from  a 
graded-index  (core)  fiber  preform,  sleeved,  and  then 
drawn  into  an  LGF.  One  purpose  of  such  a  fiber  is  to  broaden 
the  Brillouin  spectrum  in  order  to  suppress  stimulated 
Brillouin  scattering  (SBS),  among  a  number  of  other  appli¬ 
cations.’8  While  SBS  suppression  appears  to  be  a  highly 
promising  application  for  this  technology,  another  very  inter¬ 
esting  application  is  distributed  sensing  in  the  frequency 
domain.  In  a  traditional  time-domain  system,  the  Brillouin 
gain  spectrum  (BGS)  of  a  sensor  fiber  is  typically  uniform 
along  a  length  of  fiber.  Changes  in  the  value  of  vB  lead  to 
measurements  of  the  local  temperature  or  strain,  and  thus 
intuitively  it  can  be  reasoned  that  the  narrower  the  BGS, 
the  more  sensitive  the  measurement.  In  an  LGF,  however, 
since  vB  [or  more  accurately  vB(z)}  is  changing  (i.e.,  chirped) 
lengthwise  monotonically,  then  every  position  on  the  fiber 
has  a  unique  vB  with  which  it  is  associated.  If  such  a 
fiber  is  deployed  in  a  sensor  network,  then,  any  change 
locally  in  the  fiber  environment  will  result  in  a  measurable 
change  in  the  measured  BGS. 

More  specifically,  the  measured  BGS  will  be  broadened 
from  the  intrinsic  (nonLGF)  BGS,  which  is  a  Lorentzian 
function,  due  to  uB(z ),  which  also  includes  any  associated 
thermal  or  strain  distribution  along  the  fiber.  The  reverse, 
then,  is  also  true:  a  measurement  of  the  BGS  gives  rise  to 
a  measurement  of  vB(z).  If  this  distribution  is  known  for 
the  “unperturbed”  (e.g.,  room  temperature  and  zero  strain) 
fiber,  then  the  change  in  vB{z)  from  its  unperturbed  state 
gives  rise  to  a  distributed  measurement  of  temperature/strain. 
This  is  interesting  since  now  a  continuous  wave  (CW)  meas¬ 
urement  of  the  BGS,  rather  than  time-gated  returns  of  a 
pulse,  affords  a  measurement  of  the  environmental  distribu¬ 
tion  of  the  fiber,  that  is,  a  frequency-domain  system. 

In  this  paper,  the  feasibility  of  such  a  frequency-domain 
Brillouin-based  distributed  sensor  is  investigated.  The  paper 
has  three  parts.  First,  the  BGS  of  an  LGF  is  modeled,  includ¬ 
ing  the  effect  of  a  local  environmental  event.  Thermal  effects 
are  utilized  as  a  convenient  example,  where  it  is  shown  that  a 


measureable  change  in  the  BGS  is  the  result  of  a  change  in 
lengthwise  temperature  distribution.  The  corollary  case  for 
strain  or  simultaneously  both  temperature  and  strain  is 
straightforward  although  not  exemplified  herein.  Second, 
utilizing  the  LGF  optical  fiber  from  Evert  et  ah, 38  the 
model  and  the  measurable  change  in  the  BGS  are  validated 
experimentally.  Finally,  a  sensitivity  analysis  is  presented  to 
estimate  the  signal-to-noise  ratio  (SNR)  and  resolution  of 
such  a  system.  Both  pump-probe  and  coherent-detection 
schemes  for  measuring  the  BGS  are  investigated. 
Modeling  results  suggest  that  the  frequency-domain  systems 
based  on  fibers  with  a  chirped  Brillouin  frequency  have  a 
significant  potential  for  event  sensors  (such  as  alarm  sys¬ 
tems)  with  rapid  measurement  times. 


2  Modeling  the  Spectrum 

It  is  well  known  that  the  BGS  of  a  uniform  optical  fiber  is 
expressed  as  a  Lorentzian  function  of  the  frequency  v  using 
the  form 


where  vB  is  the  central  frequency  and  kvB  is  the  Brillouin 
spectral  width  (full  width  at  half  maximum).  The  central  fre¬ 
quency  is  determined  from  a  combination  of  the  Bragg  con¬ 
dition  and  the  acoustic  velocity,  Va,  via vB  =  2nmV a/X0  with 
X0  being  the  vacuum  optical  wavelength  and  nm  being  the 
effective  modal  index  as  noted  in  the  Introduction.  In 
Eq.  (1),  gBQ  is  the  peak  Brillouin  gain  coefficient,  which 
is  a  function  of  V a,  AvB,  n,  mass  density,  and  the  photoelas¬ 
tic  constant  (Pockels’  coefficient)  of  the  material  through 
which  the  light  is  propagating. 

For  LGFs,  the  central  frequency  becomes  a  function  of  dis¬ 
tance,  i.e.,  vB  ->  vB(z).  For  fibers  with  a  lengthwise  uniform 
refractive  index  difference,  this  means  that  Va  —>  Va(z)  due  to 
the  compositional  variation.  In  addition  to  the  acoustic  velocity, 
since  AvB  is  proportional  to  the  square  of  the  central  frequency 
[Eq.  (1)],  and  the  central  frequency  is  a  function  of  distance, 
Avb  is  also  a  function  of  position  along  the  fiber,  z.  Hence,  gm 
also  is  a  function  of  z  courtesy  of  Va(z)  and  A vB(z). 

The  simplest  configuration  of  the  BGS  is  a  linear  depend¬ 
ence  of  central  frequency  with  distance,  z,  along  the  fiber 
such  that  vB  — >  vB-Cz ,  where  C  is  a  constant.  In  this 
case,  the  Brillouin  spectrum  is  spread  approximately  uni¬ 
formly  across  all  available  acoustic  frequencies,  thereby  flat¬ 
tening  the  effective  BGS.  In  order  to  compute  the  measured 
BGS  from  such  a  structure,  one  assumes,  to  first  order,  that 
only  the  central  frequency  is  ^-dependent  (i.e.,  the  z-depend- 
ence  of  the  other  parameters  such  as  the  photoelastic  constant 
and  Brillouin  spectral  width  is  negligible)  and  that  the  fiber  is 
optically  lossless.  Then,  with  the  assumption  that  Eq.  (1)  rep¬ 
resents  the  local  BGS  in  the  LGF,  integration  over  all  z  (for  a 
fiber  of  length  L)  is  performed  in  order  to  obtain  for  the  mea¬ 
sured  BGS  that 


9b(v)  = 


A  vB 


2  LC 
—  arctan 


arc  tan 


,  {v  ~  vb)  +  CL 
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Fig.  1  Brillouin  gain  spectra  for  a  conventional  single  mode  fiber 
(SMF-28)  and  that  of  the  longitudinally  graded  optical  fiber  (LGF). 


along  the  fiber,  vB  will  shift  to  a  larger  value,  leaving  a  spec¬ 
tral  dip  and  a  peak  at  higher  frequency  where  the  shifted 
spectral  slice  contributes  to  spectral  slices  from  other  posi¬ 
tions  on  the  fiber. 

Modeling  the  system  constitutes  replacing  the  vB(z)  = 
vB-Cz  with  a  new,  more  appropriate  function.  For  the  exam¬ 
ple,  presented  here,  it  will  be  assumed  that  any  local  heating 
event  along  the  sensor  fiber  will  be  Gaussian  distributed  with 
some  linewidth  (spatial  extent).  Therefore,  the  new  length¬ 
wise  frequency  distribution  in  the  fiber  can  be  written  as 


vB{z)  =  vb  -  \  Cz  -  ATCt  exp 


z  -  2 

A  L 


(3) 


Next,  assumed  here  is  a  starting  central  frequency  (i.e.,  at 
one  end)  of  1 1  GHz  and  intrinsic  Brillouin  spectral  width  of 
30  MHz  to  approximately  match  the  measured  characteris¬ 
tics  of  the  SMF-28  fiber39  that  is  used  for  comparison.  The 
fiber  length,  L,  is  assumed  to  be  1000  m  with  the  additional 
assumption  that  the  fiber  has  a  linearly  distributed  total  fre¬ 
quency  breadth  of  1  GHz  (a  path  toward  such  a  fiber  using 
Ge02  and  F  co-doping  is  outlined  in  Evert  et  al.38),  then 
C  =  1  MHz/m.  It  is  noted  that  the  value  of  C  is  a  design 
parameter  that  can  be  controlled  in  the  fiber  fabrication  proc¬ 
ess.  Equation  (2)  has  been  normalized  so  that  the  total  inte¬ 
grated  Brillouin  gain  with  respect  to  v  is  conserved  between 
the  uniform  and  longitudinally  graded  fibers38,40  and  that  the 
maximum  Brillouin  gain  in  a  uniform  fiber  is  about 
2.5  X  10-11  m/W.39  Using  these  values  in  Eq.  (2),  Fig.  1 
computes  the  BGS  for  a  lengthwise  un-chirped  and  chirped 
Brillouin  frequencies.  A  clear  disadvantage  to  the  use  of  the 
LGF  in  the  Brillouin  sensing  application  is  the  considerable 
reduction  in  the  Brillouin  gain  and  the  “specialty”  nature 
introducing  a  level  of  greater  fabrication  complexity,  at 
least  currently,  over  conventional  telecom-grade  fiber. 
However,  it  is  expected  that  as  LGF  technology  matures,  fab¬ 
rication  of  such  fibers  will  become  much  more  routine. 

In  the  LGF,  a  change  in  temperature  or  strain  at  some 
position  along  the  fiber  will  result  in  a  change  in  the  local 
vB.  Usually,  as  with  Ge02-doped  silica  fibers,  an  increase 
in  strain  or  temperature  will  cause  an  increase  in  vB. 
What,  then,  happens  in  either  of  these  two  cases  is  intuitive: 
in  some  region  of  the  BGS,  corresponding  to  some  position 


where  AT  is  the  change  in  temperature,  CT  is  the  thermal 
coefficient  of  the  Brillouin  frequency  shift  (MHz/K),  and 
A L  limits  the  range  over  which  the  fiber  is  being  heated. 
Equation  (3)  is  substituted  into  Eq.  (1)  and  numerically  inte¬ 
grated  as  with  the  case  of  the  simple  linear  frequency  distri¬ 
bution.  As  an  example,  retaining  the  values  previously  stated 
and  assuming  that  A L  =  10  m,  CT  =  1.1  MHz/K  and  L'  = 
250  m  (position  of  the  heating  incident).  Fig.  2  plots  the 
BGS  for  AT  =  0,  10,  and  30  K  demonstrating  how  its 
shape  and  distribution  changes  with  a  local  heating  event. 
While  the  example  presented  here  is  that  of  heating,  similar 
results  are  obtained  for  the  strain  where  Eq.  (3)  would 
change,  so  that  AT  ->  e  (strain)  and  CT  -*■  Ce  (strain  coef¬ 
ficient  in  units  of  MHz/elongation). 

First,  changing  the  position  of  the  environmental  event 
simply  shifts  the  position  of  the  spectral  “anomaly.”  This 
is  illustrated  in  Fig.  2(a),  in  which  the  BGS  is  recalculated 
for  the  AT  =  30  K  case,  but  with  V  now  equal  to  750  m. 
Second,  since  the  intrinsic  Brillouin  spectrum  has  a  width  of 
about  30  MHz  locally  along  the  fiber,  a  reduction  in  AL 
causes  the  amplitude  or  strength  of  the  change  in  the 
BGS  to  reduce  since  a  narrower  sliver  of  spectrum  is 
being  shifted.  In  order  to  demonstrate  this,  the  calculation 
shown  in  Fig.  2(b)  is  repeated  for  AL  =  2  m,  and  the  result 
is  shown  in  Fig.  3.  Thus,  it  is  clear  that  as  the  spatial  reso¬ 
lution  (or  the  distance  over  which  the  event  takes  place)  is 
increased,  the  sensitivity  of  the  measurement  of  the  fiber 
environment  is  reduced.  A  sensitivity  analysis  is  presented 
in  a  subsequent  section. 


W) 

.S 


c 


Frequency  (GHz)  Frequency  (GHz) 


Fig.  2  (a)  Brillouin  gain  spectrum  for  the  example  for  AT  =  0,  10,  and  30  K,  (b)  expanded  view  in  the 
vicinity  of  the  thermal  activity  region  (“anomaly”).  The  amplitude  of  the  change  in  BGS  increases  with 
increasing  temperature.  Changing  L'  simply  shifts  the  location  of  the  spectral  “anomaly”  as  illustrated  with 
the  AT  =  30  K  calculation  repeated  for  750  m  [shown  in  (a)]. 
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Fig.  3  Brillouin  gain  spectrum  for  the  example  for  AT  =  0,  10,  and 
30  K  and  AL  now  equal  to  2  m. 


affording  considerable  confidence  in  the  results  and  analysis 
that  follows. 

Conceptually,  heating  and  straining  the  fiber  give  rise  to 
essentially  identical  changes  to  the  BGS.  However,  since 
heating  a  small  segment  of  the  16.7-m  fiber  was  far  easier 
to  implement  than  applying  strain,  heating  experiments 
were  performed  and  treated  here.  Heating  a  fiber  was  accom¬ 
plished  by  placing  a  segment  of  test  fiber  into  a  well-con- 
trolled  heated  water  bath.  This  type  of  thermal  distribution 
is  more  like  a  uniform  function  over  some  range  A L. 
Thus,  the  thermal  distribution  denoted  in  Eq.  (3)  is  modified 
from  possessing  a  simple  Gaussian  function  to  possessing  a 
super-Gaussian  function.  Therefore,  the  modified  final  dis¬ 
tribution  of  frequencies  for  the  LGF  is  given  by 


3  Validation  of  the  Model  and  Proof-of-Concept 

The  LGF  fiber  described  previously  in  Evert  et  al.38  is  uti¬ 
lized  in  the  following  set  of  measurements  and  modeling  cal¬ 
culations.  This  fiber,  due  to  its  acoustically  guiding  nature,41 
exhibits  multiple  acoustic  modes  in  its  BGS.  However,  the 
higher  order  acoustic  modes  will  be  neglected,  and  a 
focus  will  be  made  on  the  fundamental  (lowest  frequency) 
acoustic  mode  since  it  is  also  the  strongest.  The  measurement 
apparatus  employed  is  identical  to  that  in  Dragic,42  and  there¬ 
fore  it  will  not  be  described  in  detail  here.  However,  in  short, 
it  is  a  heterodyne  system  wherein  a  pump  signal  is  launched 
into  a  test  fiber,  generating  spontaneous  Brillouin  scattering. 
Said  signal  is  optically  amplified  and  mixed  with  the  pump 
using  a  fast  lightwave  receiver,  and  the  resulting  beat  spec¬ 
trum  is  recorded  with  an  electrical  spectrum  analyzer.  The 
measured  BGS  for  the  (F  +  Ge02  co-doped  silica  core) 
LGF  of  this  study  is  shown  in  Fig.  4.  The  spectrum  is  not 
flat,  and  thus  the  lengthwise  distribution  of  Brillouin 
frequencies  is  not  linear.  A  fit  was,  therefore,  performed 
to  determine  an  approximate  distribution.  The  distribution 
that  gave  the  best  fit  to  the  BGS  was  determined  to  be 


vb{z)=vb-{Cz+C'zXA5),  (4) 

where  vB  =  10.8985  GHz,  C  =  0.93  MHz/m,  C'  = 
1.45  MHz/m145,  and  AvB  =  39.5  MHz.  Substituting 
Eq.  (4)  into  Eq.  (1)  and  integrating  over  z  gives  the  simulated 
BGS  provided  in  Fig.  4  (as  the  dashed  line).  There  is  excel¬ 
lent  agreement  between  the  measured  and  modeled  BGS, 


Frequency  (GHz) 

Fig.  4  Brillouin  gain  spectrum  of  the  fundamental  acoustic  mode  of 
the  LGF  of  this  study  (solid  line).  The  modeling  fit  to  the  data  is 
also  shown  (dashed  line).The  total  fiber  length  is  16.7  m. 


Gi(z) 


=  vB~ 


\cz  +  C'zlAS  -  ATCt  exp 


(4r)“] 


(5) 


As  stated  before,  Eq.  (5)  is  substituted  into  Eq.  (1)  and  the 
result  is  integrated  over  z  to  yield  a  calculation  of  the  BGS  as 
a  function  of  temperature.  For  the  first  example,  L '  =  9.7  m 
and  A  L  =  1.13  m,  with  the  Brillouin  spectra  calculated  and 
measured  for  AT  =  5,  10,  20,  and  40  K.  CT  was  found  to  be 
1.11  MHz/K  for  the  LGF.5,34  Figures  5(a)  and  5(b)  show  the 
modeled  (using  M  =  20  to  simulate  the  uniform  heating  of  a 
segment  of  fiber)  and  measured  normalized  BGS  for  the  vari¬ 
ous  temperatures,  respectively.  The  normalization  procedure 
was  to  assume  the  Brillouin  gain  of  an  unchirped  (uniform) 
fiber  with  the  same  AvB  is  unity.  Figures  6(a)  and  6(b)  pro¬ 
vide  the  results  of  a  slightly  modified  set  of  conditions  where 
L'  =  9.9m  and  AL  =  1.70  m,  AT  =  5,  15,  25,  and  45  K. 
As  seen,  there  is  excellent  agreement  between  the  measured 
and  modeled  BGS.  As  described,  changing  L'  only  serves  to 
shift  the  position  of  the  dip  (“anomaly”)  in  the  spectrum; 
however,  since  the  spectral  breadth  of  the  proof-of-concept 
LGF  used  in  this  study  was  only  about  100  MHz,  performing 
such  an  investigation  would  add  little  to  the  present  analysis. 

In  summary,  it  has  been  shown  that  by  heating  (or,  in 
theory,  straining)  a  small  segment  of  LGF  (chirped  fiber) 
at  some  position,  a  measurable  change  in  the  BGS  can  be 
realized.  As  such,  in  principle,  the  reverse  process  is  then 
also  possible,  in  that  if  a  BGS  of  an  LGF  in  an  unknown 
environment  is  measured,  comparing  this  to  the  unheated 
(unstrained)  spectrum  discloses  the  thermal  (or  strain)  distri¬ 
bution.  While  the  method  used  to  measure  the  BGS  in  the 
present  analysis  is  optimized  neither  for  cost  nor  sensitivity 
and  the  available  length  of  a  fiber  was  limited,  it  provides  for 
a  reasonable  proof-of-concept.  Depending  on  the  configura¬ 
tion,  a  considerable  improvement  in  sensitivity  can  be 
achieved.  This  is  described  in  more  detail  in  Sec.  5. 


4  Limitations  to  Scheme 

4.1  Uniqueness  of  the  Brillouin  Gain  Spectrum 

For  any  known  thermal  or  strain  distribution  along  the  LGF, 
the  BGS  can  be  calculated.  This  assumes  that  the  unper¬ 
turbed  BGS  for  the  fiber  is  known.  However,  more  complex 
situations,  such  as  the  simultaneous  presence  of  multiple, 
randomly  distributed  heating  or  strain  events,  pose  a  problem 
when  a  unique  Brillouin  spectrum  from  an  LGF  is  needed  to 
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Fig.  5  (a)  Modeled  BGS  of  the  LGF  at  AT  =  0,  5,  10,  20,  40  K  and  L'  =  9.7  m  and  AL  =  1 .13  m, 
(b)  measured  BGS  for  the  same  thermal  conditions. 


identify  the  locations  of  such  multiple  events.  Although  such 
a  result  may  be  difficult  to  achieve;  in  practice,  it  is  possible 
that  two  or  more  unique  perturbation  (heating  or  strain)  dis¬ 
tributions  influence  the  BGS  in  such  a  way  that  results  in 
identical  SBS  spectra.  This  can  introduce  considerable 
uncertainty  when  the  proposed  sensing  technique  is  used 
for  applications  where  precise  identification  of  the  location 
of  an  “event”  must  be  determined.  This  will  be  illustrated 
next  through  the  presentation  of  a  simple  example. 

Consider  the  case  of  two  anomalous  heating  events  occur¬ 
ring  at  time  t  =  r0  but  located  at  different  positions  along  a 
single  LGF  (i.e.,  one  at  a  lower  and  the  other  at  a  higher  fre¬ 
quency  spectral  “position”),  as  shown  in  Fig.  7  (solid  thick 
line).  The  reader  is  reminded  that  a  spectral  position  on  the 
unperturbed  BGS  corresponds  to  a  unique  position  along  the 
LGF.  Continued  heating  for  t  >  t0  will  cause  the  spectral 
perturbations  to  shift  in  frequency  (along  the  direction  of 
the  black  arrows)  until  the  temperature  reaches  equilibrium 
or  as  heating  continues.  In  such  a  case,  the  perturbation  from 
the  heating  event  at  the  lower  frequency  position  along  the 
LGF  can  eventually  shift  in  frequency  until  it  reaches  the 
same  position  in  the  frequency  domain  as  the  second  heating 
event.  Thus,  the  spectral  slice  from  the  first  heating  event  can 
shift  and  “fill-in”  the  hole  representing  the  second  heating 
event  (as  shown  by  the  thick,  dashed  line).  The  resulting 
BGS  may,  then,  have  the  appearance  of  being  caused  by 
a  single,  much  hotter  thermal  event  with  a  larger  single 
shift  in  the  BGS,  rather  than  one  resulting  from  two  or 
more  independent  and  weaker  events.  In  such  cases,  there¬ 
fore,  multiple  heating  events  may  appear  as  fewer  heating 
events,  making  it  difficult  to  identify  the  location  and  the 
extent  of  the  various  perturbations. 


Frequency  (GHz) 


Alternatively,  the  nonuniqueness  of  the  BGS  may  intro¬ 
duce  uncertainty  in  the  actual  temperature  or  strain  distribu¬ 
tion,  even  if  the  positions  of  the  events  are  known.  This  can 
be  the  case  if  multiple  identical  events  become  spectrally 
nonsequential.  For  example,  simultaneous,  identical  events 
may  result  from  various  amounts  of  temperature  or  strain. 
In  which  case,  it  is  possible  that  some  “hotter”  or  higher 
strain  events  may  move  past  others  within  the  frequency 
domain,  making  it  difficult  to  associate  the  specific  magni¬ 
tude  of  an  event  with  its  location  along  the  LGF.  Referring  to 
Fig.  8  (for  the  example  of  fiber  heating),  the  example  BGS 
clearly  illustrates  that  at  least  two  heating  events  have  taken 
place;  however,  it  is  not  clear  which  of  the  two  is  at  the 
higher  frequency. 

4.2  Continuous  Monitoring  of  the  Brillouin  Gain 
Spectrum 

It  is  worth  noting  that  the  change  in  frequency  due  to  the  first 
event  described  in  the  first  example  above  would  need  to  be  a 
singular  value  in  order  for  the  overlapped  heating  (or  strain) 
events  to  render  the  BGS  wholly  ambiguous.  For  systems  or 
applications  where  the  possibility  of  only  a  few  anomalous 
thermal  or  strain  events  exists,  this  scenario  is  unlikely. 
However,  for  a  complex  distribution  where  perturbations 
occur  on  a  small  length  scale  and  are  very  close  together, 
this  becomes  a  more  serious  problem.  Therefore,  as  men¬ 
tioned  in  the  Introduction,  this  type  of  sensing  system  is 
best  suited  for  inexpensive  and  rapid  sense  alarms  that 
look  for  significant  “events”  rather  than  small  drifts  in  tem¬ 
perature  or  strain.  However,  in  order  to  resolve  the  problem 
of  identifying  multiple  identical  overlapping  heating  events, 
the  fiber  BGS  can  be  monitored  continuously  and  compared 
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Fig.  6  (a)  Modeled  BGS  of  the  LGF  at  A  T  =  0, 5, 1 5, 25, 45  K  and  L'  =  9.9  m  and  AL  =  1 .70  m,  (b)  mea¬ 
sured  BGS  for  the  same  thermal  conditions. 
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Fig.  7  Multiple  heating,  or  strain,  events  may  not  be  identifiable  if  the 
resulting  perturbations  to  the  SBS  spectrum  are  identical  and  overlap. 
In  this  case,  two  heating  events  appear  in  the  SBS  spectrum,  shown 
with  the  thick  solid  line.  As  temperature  increases,  the  perturbations 
on  the  SBS  spectrum  shift  with  time,  as  shown  by  the  thin  solid  lines. 
Eventually,  the  perturbations  can  overlap,  resulting  in  the  appearance 
of  a  single  large  temperature  event,  as  shown  by  the  dashed  line. 

with  those  at  previous  times.  This  enables  the  observation 
and  tracking  of  the  perturbations  within  the  SBS  spectrum 
as  they  shift  in  frequency  with  time  before  they  overlap. 
Referring  back  to  Fig.  7,  with  just  a  few  measurements  of 
the  BGS  just  after  time  t  =  f0,  but  on  a  time  scale  faster 
than  that  of  the  thermal  perturbation,  the  location  and  posi¬ 
tions  of  the  perturbation  become  clear,  despite  the  eventual 
overlapping  spectra.  Furthermore,  the  same  is  true  of  the 
example  in  Fig.  8;  measurements  of  the  BGS  as  a  function 
of  time  can  disclose  which  location  along  the  fiber  has  under¬ 
gone  a  more  significant  perturbation.  A  discussion  of  the  rap¬ 
idity,  with  which  a  measurement  can  be  made,  is  provided  in 
the  next  section. 

5  Sensitivity  Analysis 

As  previously  described,  the  system  may  take  on  two  distinct 
forms:  (1)  one  that  utilizes  a  coherent  detection  scheme 
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Fig.  8  An  example  of  a  possible  nonsequential  BGS.  It  is  clear  from 
the  spectrum  that  at  least  two  events  (local  heating  or  strain)  have 
taken  place,  but  it  is  not  readily  apparent  which  of  the  two  peaks 
(dashed  line)  corresponds  to  which  location  along  the  fiber. 
Tracking  the  BGS  in  time  can  alleviate  this  problem. 


(CDS),  and  (2)  another  that  uses  a  pump-probe  scheme 
(PPS).  The  goal  of  both  systems  is  an  accurate  and  a  very 
sensitive  determination  of  the  BGS,  from  which  can  be 
derived  the  distribution  of  a  fiber’s  environment.  Both  meth¬ 
ods  require  the  use  of  a  narrow  linewidth  spectrally  tunable 
source.  For  the  CDS,  this  is  the  local  oscillator  (LO)  and  for 
the  PPS  this  becomes  the  probe  wavelength.  The  LGF  pos¬ 
sesses  a  Brillouin  frequency  that  is  lengthwise  distributed 
and  monotonic  in  its  change.  In  principle,  this  means  that 
each  unique  frequency  component  of  the  Brillouin  spectrum 
maps  exactly  to  a  single  point  on  the  fiber.  As  such,  the  spec¬ 
tral  width  (Az/l)  of  the  laser  of  which  frequency  is  swept 
through  the  BGS,  therefore  represents  the  absolute  limiting 
spatial  resolution  for  the  frequency-domain  system  based  on 
a  LGF.  Here,  resolution  is  the  uncertainty  in  the  position 
determined  for  a  single  event.  This  limiting  spatial  resolution 
(A/)  can  be  written  in  terms  of  the  sensor  (fiber)  length,  L, 
and  total  breadth  of  frequency,  Bs,  possessed  by  the  BGS  of 
the  fiber  (for  a  linear  frequency  distribution)  as 


For  example,  assuming  a  fiber  length  of  1  km,  a  Brillouin 
spectral  breadth  of  1  GHz  along  that  fiber  length,  laser  line- 
widths  of  1  MHz  and  100  kHz  give  rise  to  absolute  minimum 
spatial  resolutions  of  1  m  and  10  cm,  respectively. 
Alternatively,  for  a  fiber  with  a  length  of  10  m  and  all 
other  parameters  kept  the  same,  the  minimum  spatial  reso¬ 
lutions  become  1  cm  and  1  mm.  respectively.  Since  the 
fiber  may  be  designed  to  be  application  specific,  LGFs 
are  not  necessarily  limited  to  long-11  or  short-range25 
applications. 

While  the  un-chirped  spectral  width  of  an  optical  fiber 
may  be  on  the  order  of  a  few  tens  of  MHz,  even  a  small 
change  in  temperature  results  in  a  measurable  and  calculable 
change  in  the  BGS.  However,  the  ability  to  measure  these 
small  changes  in  the  BGS  is  limited  by  system  noise. 
From  them,  it  is  possible  to  estimate  the  measurement  (tem¬ 
perature  or  strain)  resolution  of  a  LGF-based  system.  In  the 
following  subsections,  the  noise  sources  relevant  to  each  sys¬ 
tem  are  identified  and  listed,  and  noise  variances  (in  units  of 
Amperes  squared.  A2)  and  SNRs  are  calculated  and  pre¬ 
sented.  The  noise  variances  are  converted  into  a  measure¬ 
ment  resolution  in  the  following  way:  a  total  noise 
variance  <7t2otal  is  calculated  with  the  standard  deviation 
clearly  setting  the  minimum  detectable  change  in  the 
BGS,  as  outlined  in  Fig.  9.  The  BGS  is  calculated  for 
decreasing  temperature  or  strain,  and  the  point  where 
AA/2  (AA  is  the  magnitude  of  change  in  the  BGS  from 
the  unperturbed  stated,  as  illustrated  in  Fig.  9)  falls  to 
Gotai  is  considered  to  be  the  measurement  limit. 


5.1  Coherent  Detection  Scheme 

An  example  of  a  possible  CDS  arrangement  is  shown  in 
Fig.  10.  The  pump  signal  passes  through  an  optical  circulator 
and  into  the  sensor  fiber  through  a  second  port.  The  gener¬ 
ated  Brillouin  signal  is  backscattered  and  passes  through  the 
third  port  of  the  circulator  to  an  optical  amplifier  that  may 
possess  a  bandpass  optical  filter.  This  signal  is  then  com¬ 
bined  with  a  tunable  local  oscillator  signal  onto  a  receiver, 
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Fig.  9  Simulated  noisy  BGS  measurement  illustrating  the  standard 
deviation  of  the  noise  and  A/4.  Noises  are  assumed  to  be 
Gaussian  distributed. 


Fig.  10  Example  of  a  CDS  configuration. 


with  the  resulting  generated  electrical  signal  being  processed 
via  electronics  and  software  (not  shown). 

Two  electric  field  components  are  present  at  the  detector. 
These  are  the  LO  and  Brillouin  (designated  as  “S”)  signals: 

£lo  oc  (7) 


when  o>lo(0  =  ©s(zt).  Thus,  the  relevant  detected  homo¬ 
dyne  signal,  then,  is  2R^PLOPs. 

Next,  the  noise  variances  (units  of  A2)  are  identified  and 
listed.  The  thermal  noise  is  given  by 


4kTBN 

R 


(10) 


where  k  is  Boltzmann’s  constant  (1.38  X  10-23  J/K),  T  is  the 
absolute  temperature  (K),  R  is  a  noise  equivalent  receiver 
resistance  (Ohms),  and  BN  is  the  bandwidth  of  the  receiver 
system  (Hz).  In  the  present  system,  since  a  BGS  sweep  may 
be  1  s  or  longer,  the  thermal  and  other  noises  can  be  consid¬ 
erably  reduced  if  the  receiver  bandwidth  is  limited  to  less 
than  about  10  kHz.  Shot  noise  is  contributed  to  by  the  pres¬ 
ence  of  the  LO  and  SBS  signals,  while  the  former  usually  is 
much  larger  than  the  latter.  The  SBS  signal  is  the  total  inte¬ 
grated  SBS  power  in  the  received  signal  (over  the  entire 
BGS).  P suotai  is  distinct  from  Ps  in  that  the  latter  is  the  inte¬ 
grated  SBS  power  over  the  LO  spectrum.  Shot  noise,  then,  is 
given  by 


ffShot  —  2eR(P LO  +  P 5,total)^V  ~  2,eRP LO Bn,  (11) 

where  e  is  the  electron  charge  (1.6  X  10-19  C).  Other  dom¬ 
inant  noises  include  beat  noises  associated  with  the  various 
optical  signal  components  that  are  present.  First,  in  a  system 
that  is  optically  amplified  (as  is  the  case  with  the  CDS), 
amplified  spontaneous  emission  (ASE)  can  beat  with  the 
LO  and  SBS  signals  and  with  itself.  The  noise  variances 
are44^6 


a 


2 

ASE-ASE 


=  4N2aseB0BnR2 


(12) 


ffS-ASE  —  4(Vase(Pl o  +  ^S, total) BiyR°  »  4N  ASEPLOBNR- , 

(13) 


Es  oc  y/Fse-*”®'.  (8) 

In  Eqs.  (7)  and  (8),  P  is  the  optical  power.  The  time  or 
lengthwise  dependencies  of  the  angular  frequencies,  co, 
reflect  the  nature  of  the  signals  produced,  that  is,  the  local 
oscillator  frequency  will  be  swept,  while  the  SBS  signal 
has  a  frequency  distribution  that  is  fiber  dependent.  The 
determination  of  I^LO  +  Es\  and  the  assumption  of  a  detec¬ 
tor  responsivity,  R  (A/W),  gives  rise  to  the  following  current 
signal  that  is  produced  by  the  system 

4ig  =  r(Ax>  +  ^s  +  2  \/PloPs  cos{  V°lo  (0  —  ■ 

(9) 

In  Eq.  (9),  although  Ps  may  not  be  polarized  identically 
with  PLO,  it  is  assumed  in  the  subsequent  analysis  that  the 
polarization  components  can  be  treated  independently  (opti¬ 
cally  polarization  split)  and  that  the  resulting  electrical  sig¬ 
nals  are  added  after  the  electrical  mixing.43  This  is 
considered  to  be  a  reasonable  first-order  approximation 
given  in  the  exploratory  nature  of  this  sensitivity  analysis. 
In  a  homodyne  system,  a  change  in  the  direct  current 
level  is  detected  as  coLo(0  is  swept  through  cwSBS(z)  or 


where  B0(»  BN)  is  the  optical  bandwidth  of  the  ASE  noise 
spectrum,  /VA se  (W /Hz),  which  is  assumed  to  be  optically 
filtered  with  a  passband  of  B0.  Similarly,  since  there  is  also 
Brillouin  spontaneous  scattering,  these  beat  noises  may  be 
presented  as  well.  Similar  to  Eqs.  (12)  and  (13),  including 
an  ASE-SBS  term,  these  can  be  described  as 


n  2 

°SBS-SBS 

=  4N^B  SB  NR2 , 

(14) 

°LO— SBS  " 

=  4  NSP  loBnR~  . 

(15) 

_2 

^ASE-SBS 

=  4NSN  a$eBsBnR- , 

(16) 

where  Ns  is  the  spectral  power  density  (W/Hz)  of  the 
generated  Brillouin  signal,  and  Bs  is  the  spectral  breadth 
of  the  fiber  BGS  (units  of  Hz). 

For  a  uniformly  distributed  BGS,  the  Brillouin  signal 
power  is  related  to  the  power  density  and  LO  laser  linewidth 
as  Ps  =  NsI\vl.  In  Eq.  (16),  it  is  assumed  that,  most  typi¬ 
cally,  Bs  <S  B0.  Seven  noise  sources  have  been  identified, 
and  other  possible  noise  sources  (such  as  laser  relative  inten¬ 
sity  noise)  can  be  rendered  negligible  through  careful 
selection  of  system  components. 
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In  order  to  calculate  Ps(=  AvsNs),  the  analysis  found  in 
Boyd  et  al.47  was  utilized  where  a  Brillouin  reflectivity,  RB, 
was  used  to  determine  the  total  Brillouin  signal  density, 
NS(NS  =  PRb/Bs)  appearing  at  the  receiver, 


Rb  =  YeG'l 2 


(17) 


where  is  the  modified  Bessel  function  of  order  m  and 


Y  ={h+  \)gBhvsr—, 

(18) 

G'  =  gBjL, 

(19) 

(hv, A  ]-> 

(20) 

and  where  P  is  the  pump  power  launched  into  the  fiber,  A  is 
the  mode  effective  area  (determined  from  the  MFD),  vs  is  the 
frequency  of  the  Stokes’  signal  (approximately  equal  to  that 
of  the  pump  signal),  and  F  is  the  decay  rate  (r  =  nBs). The 
Brillouin  frequency  vB  was  arbitrarily  selected  for  conven¬ 
ience  to  be  11  GHz,  and  the  remaining  assumed  values 
can  be  found  in  Table  1 .  Since  Ns  is  small,  it  is  anticipated 
that  the  optical  preamplification  will  be  applied  via  an  EDFA 
and  an  optical  power  gain,  G,  of  30  dB  is  assumed  available. 
Thus,  in  each  of  the  relevant  noise  equations,  the  fiber-ampli¬ 
fied  value  GNS  is  utilized.  Finally,  the  total  SNR  of  the 
system  can  be  calculated  using 


SNR  = 


2R^7W^ 

v/EM 


(21) 


where  the  summation  is  over  the  identified  noise  variances 
[Eqs.  (10)— (16)]. 

Next,  an  example  is  presented  to  determine  the  relative 
strengths  of  each  of  these  sources  of  noise.  Table  1  lists 
the  assumptions  made  for  the  example  of  a  distributed  tem¬ 
perature  sensor.  The  optical  bandwidth  Ba  is  assumed  to  be 
conditioned  by  a  bandpass  filter,  such  as  a  Bragg  grating 
assembly,  and  the  SBS  bandwidth  Bs  is  assumed  to  be  lim¬ 
ited  by  the  spectral  extent  of  the  fiber  BGS,  in  this  case, 
assumed  to  be  1  GHz  linearly  distributed  over  a  fiber  length, 
L ,  of  1000  m,  such  that  the  BGS  is  approximately  uniform 
in  shape. 

Figures  11  and  12  display  the  logarithm  of  the  calculated 
noise  variance  for  each  identified  source  (for  BN  =  1  kHz) 
and  the  SNRs  versus  P  for  several  BN,  respectively.  Both  the 
thermal  and  shot  noises  are  dominated  by  the  other  noises 
and  can  be  neglected  here.  Inspection  of  Fig.  12  shows 
that  there  is  a  range  of  optimal  P' s  where  the  SNR  is  maxi¬ 
mized.  This  can  be  explained  as  follows.  As  the  power 
increases,  SNR  increases  as  expected.  However,  since  the 
second  dominant  noise  source  [SBS-SBS  Eq.  (14)]  is  depen¬ 
dent  on  N2s  and  Ps  on  N$,  at  very  large  values  of  P,  the  SNR 
will  begin  to  decrease  as  the  inverse  of  NlJ2. 

In  order  to  estimate  the  sensitivity  or  thermal  resolution  in 
the  context  of  the  previous  set  of  calculations,  the  method 
outlined  in  the  beginning  of  this  section  is  utilized.  Strain 


Table  1  Values  assumed  for  the  example  calculation  presented 
here. 


Parameter 

Symbol 

Value  assumed 

Unit 

Local  oscillator  power 

Plo 

1 

mW 

Detector  responsivity 

R 

1 

A/W 

System  electrical  bandwidth 

Bn 

0.1,  1,  10 

kHz 

System  temperature 

T 

300 

K 

Noise  equivalent  resistance 

R 

1 

kn 

Center  optical  wavelength 

Aq 

1550 

nm 

Optical  passbandwidth 

Bo 

0.1 

nm 

Optical  passbandwidth 

B0 

12.5 

GHz 

Total  Brillouin  spectral  width 
of  the  LGF 

Bs 

1 

GHz 

Sensor  length 

L 

1000 

M 

LO  or  probe  laser  spectral 
width 

Ai jl 

1 

MHz 

Brillouin  frequency  shift 

VB 

11 

GHz 

Mode  field  diameter 

MFD 

10 

/yiti 

Fiber  amplifier  gain 

G 

30 

dB 

Spontaneous  emission  noise 
power  spectral  density 

N/kse 

io-15 

W/Hz 

Brillouin  gain  coefficient 

9b 

0.12x10-” 

m/W 

Brillouin  shift  temperature 
coefficient 

Ct 

1.11 

MHz/K 

Brillouin  shift  strain  coefficient 

c, 

500 

MHz/% 

is  considered  in  the  next  paragraph.  Assuming  the  resolution 
limit  for  the  present  case  (AvlL/Bs  =  1  m)  and  an  SNR  of 
around  100  (100  Hz  BN  case),  a  temperature  resolution  of 
about  10  K  is  obtained.  While  at  first  glance  this  may  not 
seem  very  impressive,  this  system  resolution  is  for  a  single 
BGS  sweep  that  can  be  performed  as  quickly  as  10  ms  due  to 
Bn  (and  limited  by  the  electronics  and  laser  tuning  speed). 
As  such,  if  a  measurement  integration  can  be  done  such  that 
the  measurement  of  a  thermal  distribution  is  averaged  over 
several  seconds,  the  noises  can  be  greatly  diminished,  and  a 


Fig.  11  Logarithm  of  the  variances  of  the  various  identified  noise 
sources. 
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Fig.  12  SNR  of  the  homodyne  system  plotted  as  a  function  of  P  for 
three  different  electrical  bandwidths  BN. 

considerable  improvement  in  the  sensitivity  can  be  achieved. 
This,  of  course,  will  depend  on  system  requirements,  and 
therefore  a  detailed  analysis  is  not  provided  here.  As  it 
stands,  a  single-sweep  system  could  be  implemented  as  a 
rapid-detection  thermal  alarm  system  for  a  structure. 

Finally,  this  subsection  is  concluded  with  a  discussion  of 
the  strain.  In  principle,  the  measurement  of  strain  rather  than 
temperature  is  entirely  feasible,  and  the  entirety  of  the  above 
noise  analysis  holds  for  that  case.  However,  in  the  sensitivity 
analysis,  the  assumption  is  made  that  the  Brillouin  frequency 
shift  increases  at  a  rate  of  about  500  MHz/%  of  strain.  Thus, 
recalculating  the  strain  sensitivity  using  the  identical  set  of 
assumptions  provided  above,  a  minimum  detectable  applied 
strain  is  found  to  be  about  0.02%.  Again,  this  is  the  sensi¬ 
tivity  in  the  lengthwise  resolution  limit  in  a  single  measure¬ 
ment  sweep. 

5.2  Pump-Probe  Scheme 

An  example  of  a  possible  PPS  arrangement  is  shown  in 
Fig.  13.  The  pump  signal  passes  through  an  optical  circulator 
and  into  the  sensor  fiber  through  a  second  port.  A  probe  sig¬ 
nal  is  launched  into  the  far  end  of  the  sensor  fiber  and  expe¬ 
riences  gain  due  to  SBS  in  the  sensor  fiber.  The  Brillouin- 
amplified  probe  signal  then  passes  through  the  third  port  of 
the  circulator  and  into  a  receiver.  The  resulting  generated 
electrical  signal  is  processed  via  electronics  and  software 
(not  shown). 

In  the  case  of  the  PPS,  the  measured  signal  is  that  of  an 
amplified  probe  source  that  backpropagates  through  the  sen¬ 
sor  fiber  in  a  direction  opposite  to  that  of  the  pump.  Gain 
originates  from  SBS  itself,  and  thus  the  need  for  optical 
amplification  generally  is  deemed  unnecessary.  The  ampli¬ 
fied  probe  power  at  the  receiver  can  be  expressed  as  a 


Fig.  13  Example  PPS  configuration. 


function  of  the  Brillouin  gain  provided  in  Eq.  (19)  and 
the  input  probe  power  (Pprobe)  as  P's  =  PpmbeeG  ■  From  a 
system  perspective,  an  unamplified  probe  signal  (outside 
of  the  Brillouin  gain  band)  leads  to  an  offset  in  the  measured 
signal,  much  like  a  nonzero  space  (a  zero  in  a  binary  system) 
in  a  lightwave  communications  system.  As  a  result,  mean¬ 
ingful  information  in  the  measured  Brillouin-amplified  sig¬ 
nal  is  therefore  the  change  in  power,  which  gives  rise  to  an 
electrical  signal  (current)  having  the  form 

*Sig  =  R^probe(eG'  -  I)-  (2D 

Since  an  optical  amplifier  is  not  used,  all  but  Eqs.  (12), 
(13),  and  (16)  remain  relevant  to  the  PPS  configuration.  In 
moving  to  the  PPS  from  the  CDS,  Eqs.  (10)  and  (14)  remain 
the  same,  while  the  optical  power  terms  in  Eqs.  (1 1)  and  (15) 
are  replaced  by  P's  [i.e.,  PL0  +  P$  P's  and  P|.o  P's  in 
Eqs.  (11)  and  (15),  respectively]  since  P's  represents  the  total 
optical  power  at  the  receiver.  The  assumption  is  made  that  the 
Brillouin  spontaneous  scattering  signal  is  unchanged  with 
changes  in  the  probe  power  (i.e.,  un-depleted  pump  approxi¬ 
mation),  which  is  reasonable  at  low  probe  power. 

Utilizing  a  pump  power  of  100  mW,  Figs.  1 1  and  12  are 
repeated  for  the  PPS  configuration  as  a  function  of  probe 
power  and  are  provided  as  Figs.  14  and  15  utilizing 
BN  =  1  kHz  and  identically  all  relevant  quantities  found 
in  Table  1.  From  the  figures,  it  is  clear  that,  at  least  in 
terms  of  the  SNR  (=  is ig/ \/X^=  i  the  PPS  is  vastly 

superior  to  the  CDS.  This  can  be  understood  by  considering 
that  (1)  the  signal  power  in  the  PPS  scheme  is  much  larger 
than  in  the  CDS  scheme,  (2)  the  spontaneous  Brillouin  scat¬ 
tering  noise  is  not  preamplified  before  entering  the  receiver, 
and  (3)  the  gain  term  in  Eq.  (22)  prevents  the  Brillouin  spon¬ 
taneous  scattering  noise  from  dominating  the  system  at  any 
pump  power.  The  advantage  of  the  CDS,  however,  over  the 
PPS  is  that  it  is  a  true  reflectometry  system,  whereas  the  PPS 
requires  either  a  probe  laser  at  the  far  end  of  the  sensor  fiber 
or  some  type  of  reflector  (such  as  a  coated  end-facet  or  fiber 
Bragg  grating)  at  the  far  end  of  the  sensor  fiber  or  a  second 
fiber  and  loop  through  which  one  launches  the  probe  signal. 
Any  of  these  possibilities  may  be  infeasible  depending  on  the 
system  or  environment. 

The  SNRs  utilizing  other  electrical  filter  bandwidths  were 
not  calculated  due  to  the  already  large  values.  Using  the 
results  in  Fig.  15,  and  assuming  an  SNR  above  103  can 
be  achieved,  the  sensitivity  of  the  system  can  be  estimated 


Pprobe  (W) 


Fig.  14  Logarithm  of  the  variances  of  the  various  identified  noise 
sources  for  the  PPS  configuration.  The  variance  is  calculated  in 
units  of  Amperes  squared. 
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Fig.  15  SNR  of  the  pump-probe  system  plotted  as  a  function  of  probe 
power  for  BN  =  1  kHz. 

utilizing  the  same  method  as  above,  giving  rise  to  one  that  is 
below  0.5  K  in  a  single  frequency  sweep  of  the  probe  laser. 
The  lengthwise  absolute  resolution  of  the  system  in  this  case 
remains  1  m,  as  with  the  CDS  example,  limited  by  the  probe 
laser  spectral  width.  This  suggests  a  highly  sensitive  and 
rapid  measurement  system  is  possible  in  the  frequency 
domain  utilizing  LGF  in  a  chirped  fiber  sensor  configuration. 

6  Conclusion 

A  new  type  of  frequency  domain,  Brillouin-based  distributed 
sensing  system  was  proposed,  outlined,  and  evaluated.  The 
system  utilizes  a  LGF,  of  which  Brillouin  frequency  changes 
with  position  along  the  fiber,  which  then  is  an  intrinsically 
“chirped  fiber.”  Since  each  frequency  on  the  Brillouin  spec¬ 
trum  maps  one-to-one  to  a  position  along  the  fiber,  a  change 
in  the  Brillouin  spectrum  at  some  frequency  represents  an 
environmental  event  at  the  same  position  along  the  fiber. 
As  such,  the  distributions  of  temperature  or  strain  will  result 
in  a  measureable  change  in  the  whole  BGS.  Experimental 
measurements  on  an  LGF  with  roughly  100  MHz  acoustic 
frequency  chirped  distributed  over  16.7  m  demonstrated 
the  potential  of  such  a  frequency-domain  system. 
Uniformly  heated  were  1.1  and  1.7  m  segments  of  this 
LGF,  through  which  a  model  was  shown  to  extremely 
well  predict  the  measured  BGS,  rendering  the  BGS 
calculable  for  a  known  temperature/strain  distribution.  The 
segments  were  uniformly  heated  up  to  40  K  above  ambient. 

In  a  system,  with  knowledge  of  the  unperturbed  spectrum,  the 
perturbed  fiber  spectrum  gives  rise  to  a  lengthwise  measurement 
of  the  temperature  or  strain.  Since  the  system  requires  CW  lasers 
and  a  BGS  is  measured  in  order  to  determine  T(z)  and  e(z),  this 
can  be  classified  as  a  frequency-domain  distributed  sensing 
system.  Sensitivity  analyses  based  on  coherent  detection  and 
PPS  were  performed  and  showed  that  these  systems  have  the 
potential  to  achieve  very  high  SNR  values  with  sub-1  K  and  sub¬ 
meter  measurement  accuracy  in  a  single,  extremely  rapid  scan  of 
the  BGS.  Pump-probe  methods  seem  to  have  a  considerable  ad¬ 
vantage  over  their  coherent  counterparts. 

The  mathematical  analysis  presented  here  was  limited  to 
the  case  of  distributed  sensing  of  events  that  may  trigger  an 
alarm,  such  as  a  fire,  intrusion,  or  leak.  High  spatial  resolu¬ 
tion  with  the  ability  to  resolve  two  closely  spaced  simulta¬ 
neous  events  will  require  LGFs  with  greater  total  frequency 
span,  or  larger  values  of  C  (MHz/m).  Furthermore,  measure¬ 
ments  were  limited  due  to  the  relatively  short  available 
length  of  proof-of-concept  fiber.  As  such,  work  is  currently 


underway  to  increase  the  length  of  fibers  and  expanding  the 
total  breadth  of  the  Brillouin  spectrum,  and  results  are 
forthcoming. 
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